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The most difficult problem in shallow underwater acoustic communications 
is considered to be the time-varying multipath propagation because it impacts 
negatively on data rates. Computationally intensive and complex signal 
processing algorithms are required to compensate for symbol overlapping. This 
thesis presents results of a tank scale experiment to test Time-Reversal 
Acoustics (TRA) approach for high data rate binary transmissions. TRA can 
environmentally adapt the acoustic propagation effects of a complex medium. 
Our results show the suitability of the TRA approach in underwater acoustic 
communications. The results also show good focusing properties at an intended 
target location. The focal region extends over a few wavelengths, outside of 
which scrambling of the message occurs, offering natural encryption. Range 
shifts of the focal region could be achieved by frequency shifting. We found that 
the time focusing is aperture-size independent, but the spatial focusing is 
aperture-size dependent. Overall, we showed that our algorithm can accomplish 
a fast, secure, and stable communication scheme with low computational 
complexity. 
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I. INTRODUCTION 



For underwater acoustic propagation in shallow water environments, a 
narrow-band pulse generated by a point-like source spreads because of multi- 
path propagation due to reflections at the surface and at the bottom of the ocean 
or refractions within the interior. Each propagation path takes a different time, so 
that the transmitted pulse gives rise to multiple copies, not equally spaced in 
time, arriving at a receiver (Fig. 1.1). The strong signal degradation due to the 
multipath propagation is further enhanced by high spatial and temporal variability 
of the channel conditions that exist in a typical underwater acoustic channel in 
shallow water environments. 





transmitted signal 



received signal 



Figure 1.1. Example of transmitted and received signal in the multipath propagation 
environment of a rectangular waveguide. The time-scale for both signals is the same. 



In a communication link, the signal degradation due to multipath spreading 
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causes Inter-symbol Interference (ISI) which is the interference between symbols 
due to multiple arrivals of each transmitted symbol. For example, in a shallow 
water channel a multipath spread of 100 ms with 1000 symbols per second, the 
ISI extends to over 100 symbols. Also, as long propagation times decrease data 
throughput, available bandwidth is severely limited by transmission loss, which 
increases both in frequency and in range. 

A novel technique called Time-Reversal Acoustics (TRA) can 
environmentally adapt the acoustic propagation effects of a complex medium in 
order to focus energy at a particular target range and depth. Additionally, 
estimates show that TRA focusing can be done at useful ranges (30 km) and 
over reasonable time scales (30 min) without updating the time-reversed signal 
under certain conditions, as shown by Kuperman et al., (1998) in an ocean 
experiment using 445 Hz frequency signals. The TRA technique consists of 
digitizing the analog signal received by, say, a microphone or a hydrophone, time 
reversing it, and retransmitting it from a co-located source. If a wide-aperture 
array of receivers/transmitters is used, the time-reversed signal back at the 
point-like source is focused in time and space. For a small-aperture array, while 
the time-reversed signal back at the point-like source does not focus well In 
space, it still exhibits near-optimal focusing in time. 

Because the path back from the array to the original source location is 
reciprocally identical (barring significant currents or other features which degrade 
reciprocity), the mode group speeds are also identical. Therefore, all of the 
modes arrive back at the original source location at the same time and add up 
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constructively. This generates both a temporal and spatial (range and depth) 
focusing of the propagating field. Thus, using TRA the multipath structure is 
eliminated because all of the propagation paths add coherently at the source 
location, which results in a reduced ISI of the communication link. Furthermore, 
due to spatial and temporal focusing there is a strong enhancement of the field at 
the focus location, which allows us to increase bandwidth and the transmission 
range. Also, because the spatial extent of the focus is a few wavelengths of 
sound, TRA possesses natural encryption for points other than the intended 
target. 

An important step toward useful naval and commercial applications of 
TRA is to fully explore the capabilities and limits. This exploration, often 
repetitive in nature, is most cost effectively done in tank-scale experiments, and 
supplemented with numerical simulations. 

The properties of time-reversal acoustics suggest a potential application in 
the field of non-coherent acoustic communications. This has previously been 
suggested by Jackson and Dawling (1991) and Kuperman et al. (1998). “Non- 
coherent” in this sense is meant to imply that no phase information will be 
transmitted in the communication signal but instead an energy detector method 
will be used. A novel communications scheme has been suggested by Abrantes 
et al. (1999) and Smith et al. (1999), which employs both the temporal and 
spatial focusing properties of TRA. In this scheme, two source/receiver systems, 
A and B, are used to communicate with one another. A set of distinct signals 
representing binary bits is agreed upon and known by the A and B systems. 
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Each signal should have roughly the same temporal resolution. A simple 
example would be multiple broadband signals with identical bandwidth, Hanning 
window source spectra, and center frequencies separated by half the bandwidth 
(allowing for null detection between adjacent signals). The total number of 
signals within the set is then limited primarily by the bandwidth of the system. 
For this example, assume that there are four distinct signals (bits) within the set 
defining a 4-bit character (giving values from 0-15). The communications link is 
achieved when a single element of A transmits each signal to B with sufficient 
time between transmissions for all multipath structure to arrive. This provides B 
with the transfer function of the environment for those signals. Then B has the 
ability to transmit binary information to A by building 4-bit “symbols” out of the 
time-reversed signals. In addition, the transmission of subsequent symbols can 
overlap each other significantly. The criterion for the delay between subsequent 
symbol transmissions is the temporal resolution of the focus back at A. For 
example, if the signals originally transmitted by A have a nominal temporal 
resolution of 1 0 msec, then B is able to transmit the 4-bit symbols at a nominal 
rate of 100 symbols/sec. Furthermore, recall that these focused signals arriving 
back at A have an improved SNR. Now consider that multiple elements in A 
each uniquely transmit the communication signals. Depending upon the spatial 
focusing properties of the time-reversed signals, multiple symbols may be 
transmitted from B to focus at different locations along A. Thus for a full 
transmission bandwidth of 1 kHz and a pulse bandwidth of 200 Hz, this scheme 
can provide a rate of 2,000 bits/s. This is comparable to current rates with the 
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added benefits of longer range due to energy focusing and covert coding due to 
the inherent scrambling induced by the environment at points other than the 
intended receiver. 

The aim of this thesis is to establish an experimental study for the results 
found in previous simulations and to provide further knowledge about the relation 
of window size versus focal region size. Also the practicality of the method in 
underwater acoustic communications will be shown. 

The organization of this thesis is as follows. In Chapter II, we describe the 
theory for TRA. Chapter III deals with the description of the apparatus. Chapter 
IV gives a description of the experiment and the definition of the symbols. 
Chapter V presents the results of different TRA measurements and the 
conclusions. 
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II. TIME REVERSAL ACOUSTICS THEORY 



This chapter presents an overview of time-reversal acoustics (TRA) 
theory. We concentrate here on the pulse excitation because of the relevance it 
bears with our experiment. Also we discuss the focal region properties and the 
issues of moving communications partners. 



A. PULSE EXCITATION 

In the following we give an overview of the theory for pulse excitation as 
described in Kuperman et al. (1998). For the geometry as shown in Fig. 2.1 we 
have for the element of the receiving array, consisting of a total of J. elements, 
the following time-domain signal 

P(i?,z^.,r) = J GJR,Zj,z,)S(oi)e-‘‘“d(o (2.1) 

Here P(R,Zj,t) is the pressure field as seen by the array element at 
range R, depth Zj and time t. G^(R,Zj,z,) denotes the transfer function of the 

environment at the position of the element of the receiving array due to a point 
a point source at range r=0 and depth z=Zs and S(ty) is the initial pulse 
transmitted at frequency co. This expression includes all effects of the waveguide 
environment, such as time elongation of the pulse due to multipath propagation. 
For convenience we set the time origin such that P(R,Zj,t) = 0 outside the time- 

interval [0,t] , where t is a time greater than the length of the time-stretched 
signal. The time-reversed signal that will be retransmitted from the j“’ element will 
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then be P{R,Zj,T-t), where r>2ris an arbitrary point in time. This condition is 

imposed by causality. The signal has to be completely received before it can be 
time reversed. Note that the time reversal in the time-domain is the equivalent to 
the phase conjugation in the frequency-domain. The time-reversed pressure field 
can be written as 



Thereby the sign of the integration variable co has been reversed and the 
conjugate symmetry of the frequency-domain Green’s function and initial pulse 
has been used. * denotes the complex conjugate. The bracketed expression in 
Eq. (2.2) is the Fourier transform of the signal received by the j*^ array element 
after time reversal and time-delay. Further we find that this expression is the 
frequency-domain representation of the retransmitted signal. We can therefore 
use Fourier synthesis to obtain the time-domain representation of the field 
produced. The resulting acoustic field satisfies the wave equation 



where k{z) is the depth dependent wave number. 

For a vertical line of discrete sources we find the solution of Eq. (2.3) - 
using Green’s function theory - to be 






( 2 . 2 ) 




J 



V*/^,(r,z) + r*(z)/;,(r,z) = X<^(z-Z;)G:(/?.Z;.z,) 



(2.3) 




(2.4) 
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Equation (2.4) can be used to show that the time reversal process 
produces focusing in time as well as in space. The focusing in time occurs 
because a form of matched filtering is done. 

To show this, we examine the time revered field at the focus point, 
meaning we set in Eq. (2.4) r = R, z = z, . If we neglect the density gradients, 
reciprocity allows the interchange G^(R,z,,Zj) = G^(R,Zj,z,) (Jensen et al. 

(1993)). Using the time-domain representation of the Green’s function and initial 
pulse the time-domain equivalent of Eq. (2.4) becomes (Kuperman et al. (1998)) 



P(r,z,r) = — G,.„.(/?,z,.,z,)xG,.(/?,z,,z,)Jr’ls(r"-r + r)Jr" (2.5) 
{27t) j=i L'' J 



Note that the Green’s function is correlated with itself, which is a matched 
filtering operation. The time elongation due to multipath propagation is reduced, 
hence we focus in time. The sum over all elements form a spatial matched 
filtering, that is to say we focus also in space. 



9 



R 




Figure 2.1. The geometry of the experiment as it is used in the description of the time 
reversal theory. 
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B. FOCAL REGION PROPERTIES 



The primary factors influencing the time reversal focusing are the size of 
the time window chosen and the degradation of the signal due to the 
environment. The later can be divided into static and dynamic categories. The 
static category includes propagation and array structure effects. The dynamic 
one includes the time-varying environmental effects (Kuperman et al., 1998). 

We can write the field at the array in the general form 

/^,(r,z,ty) = XG3(r,z;0,zpG,*(0,z^;i?,zJ (2.6) 

Thereby G, (0,Zy;i?,Zj) is the Green’s function describing the propagation 
form source to array and G 2 (r,z; 0 ,Zy) the one describing the propagation from 

array to any field point. The subscripts 1 and 2 allow for a changing environment 
within the two transmissions. However, we assume that during each transmission 
the environment is "frozen". 

1. Time reversal focusing in static environment 

We will examine the time reversal focusing in static environments by 
considering a general nonuniform, nonadiabatic waveguide. First we will derive 
the "ideal" case and then identify the degrading factors. For simplification we 
consider only vertical arrays. We assume small regions to be range-independent 
but larger regions to be range-dependent. 
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The Green’s function is approximated using range-independent normal 
modes. The Green’s function near the source can therefore be written as 
(Kuperman et al., (1998)) 



G^{r,z\R,z,) = 

r^k„(R)\r-F\ 



(2.7) 



Similarly the Green’s function near the array is 

n (c\ R ^ .<=.( 0 )/? 



( 2 . 8 ) 



The modal eigenfunctions are denoted as u^(R,z) and u„(0,Zj), 
respectively. The corresponding eigenvalues are k„(R) and /:„(0). These 
Green’s functions do not carry the subscripts 1 and 2, because we here look at 
the time-invariant case. The mode amplitudes for the near source and near array 
expressions are (Kuperman et al., (1998)) 



and 










(2.9a) 



Kiz) = 'Z^„^a„(z). (2.9b) 

n 

p is the density of the medium. 

For convenience we assume that there are the same number of modes 
excited at both locations. This leads to being a square matrix (Kuperman et 
al., (1998)). The matrix includes any mode coupling due to the range- 
dependence and is defined independent of the source depth. Also Is unitary. 
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as far as the absorption loss in the water column and into the bottom can be 
neglected. 

The field produced by the array is described by Eq. (2.6), if we do not 
allow for a time-variant environment. The Green’s function for the propagation 
from the element to any field point therefore, in terms of mode amplitudes, 
becomes 



G^{r,z\0,Zj) 



P{Z) y 



( 2 . 10 ) 



Where the mode amplitudes c„(z) are essentially the same as b„(z) , but with the 
source range shifted by r-/? 



c„(z) = ^U^a„(z)e“‘’ 



(RXr-R) 



( 2 . 11 ) 



Putting everything together we get an expression for the time-reversed field in a 
range-dependent waveguide (Kuperman et al., (1998)) 



P,^{r,z,(o) 



Pjz) y c„,(z)A„„,^?;(z,)^. /K(0)-t;(0)1/? 
^|kJ0)k:(0) 



( 2 . 12 ) 



Here is 

nvi 



^ y u„(0,Zj)u„(0,Zj) 

Piz,) 



(2.13) 



In the ideal case, the array spans the full water column and its elements 
have equal spacing . In the bottom the amplitudes are negligible. In this case 
the sum over array elements in (2.13) approximately equals the orthogonality 
integral for modal eigenfunctions. For this case we get (Kuperman et al., (1998)) 
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(2.14) 



P,,{r,z,(o) 



Pi.z) 

d.4Tr 



m,n 



where 



I 



^Im^ln .-ZlMk.mR 

KiO) 



(2.15) 



The losses due to scattering and absorption are detrimental to time- 
reversal focusing because they attenuate higher-order modes, resulting in a 
blurrier focus than without loss. This blurring will increase with distance between 
source and array due to the strong range and mode number dependence of 
attenuation. 

For the ideal case, that is no losses and unitary mode coupling matrix , 
we find 



(2™,= 



Km 



(2.16) 



This leads to a time-reversal field for an ideal array in a lossless environment of 
the form (Kuperman et al., (1998)) 






r SnpiZs)k„iO)d,yf^ 



(2.17) 



Even though Eq. (2.17) represents the ideal case, it illustrates the 
properties of actual time-reversal arrays, as long as they are not too far from 
ideal. The independence of the focus region from the distance between source 
and array, as long as losses can be neglected, is one such property. This means 
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that only local environmental properties affect the focus region and that there is 
no range dependence. 

The derivation above is essentially two-dimensional. If static out-of-plane 
scattering occurs, this will lead to degradation if one-dimensional arrays are 
used. However, this will not affect planar or spherical arrays. 

2. Time reversal focusing in time-varying environments 

In a real environment there occur changes in time due to surface and 
internal waves and sound speed profile changes. To discuss their influence on 
time reversal focusing it is convenient to decompose the Green’s function into a 
coherent and an incoherent part (Kuperman et al., 1998), 



In Eq. (2.18), a is either 1 or 2, the subscripts as used before. The coherent part 
G is without subscript because it is unchanging. We assume that there is 
an uncorrelated change in the environment. 



Where 7^,7^ and ^ are the vectors to arbitrary points in space. Using Eqs. (2.6), 
(2.18) and (2,19), we find the mean time-reversed field to be 



G„(F,r’) = G(F,F’) + (5G„(r,r’) . 



(2.18) 




(2.19) 



J 



P,^{r,z,(0) = Y^G (r, z;0, z^G * (0, Zj\R,Z,) 



( 2 . 20 ) 
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where " denotes the averaged value and the subscript tr stands for time- 
reversed, and the variance of the field 



P,^{r,z,(of -\P„{r,z,(of 






G (r, z;0, Zj)G’ (0, Zj ;R,z,)Kjj,(R,Zs)— 

+G iR,z/,0,Zj)G\R,z,;0,Zy)Kjj.(r,z ) + K^.{r,z )Kjj.{R,z,) 



( 2 . 21 ) 



where 



K^.{r,z) = {SG,{0,zf,r,z)8Gl{Q,Zy,r,z)) . (2.22) 

The covariance Kj^.(j,z) is proportional to the correlation between the incoherent 

field at the j*'’ and element of the array. We used reciprocity properties to get 
Eq. (2.20), allowing for interchange of the Green's function arguments, and 
stationarity, meaning that <5G, and <5^2 have the same statistics. 

Note that Eqs. (2.20) and (2.21) are general and include three- 
dimensional cases. They lead to the conclusion that the mean focus field can be 
obtained by using the coherent Green's function rather than the actual, random 
one. Further it is found that the field near the source is not fluctuating 
significantly. 
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3. Time reversal focusing range shifts 

Song et al. (1998) stated that the sound field in an ocean environment is 
characterized by the existence of lines of maximum intensity I(w,r) having a 
slope This slope is invariant for a particular group of modes. 



Here v and u are the phase and group velocities characterizing the given mode 
group. It is further stated that is an invariant because it is determined by the 
properties of the medium and the derivative in Eq. (2.23) averaged over the 
mode group of interest. The invariant property applies to the sound intensity, the 
phase velocity and the group delay. Therefore the invariant property holds for 
broadband pulses or multiple center frequency narrowband signals, provided that 
the mode shapes do not change significantly over the bandwidth. 

Song et al. (1998) derived that a frequency shift leads to a shift of the focal 
range. The new focal range K can be expressed relative to the old focal range 
as 



For the case of a Pekeris waveguide with this leads to a relation where an 
increase of the frequency by p percent leads directly to an increase in range of p 
percent. 




(2.23) 




(2.24) 
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It should be noted that Eq. (2.24) is an approximation in the sense that the 
frequency dependence of the mode amplitude is neglected. Also the expansion is 
only valid for a bounded group of modes for a given R and co, and further 
multimode propagation was assumed in the derivation. 
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III. APPARATUS 



In this chapter we describe the apparatus, instrumentation, and 
experimental techniques used to communicate in a shallow water environment 
using acoustic time-reversal techniques. The main components of the experiment 
were the computerized data generation and acquisition system with its two 
transducer arrays. 



A. THE ACOUSTIC ENVIRONMENT 

The experiment is conducted in a fiberglass coated wooden tank. The 
side-walls of the tank are covered by anechoic tiles to remove the influence of 
cross-range reflections. This tank simulates a shallow-water environment 
waveguide with infinite horizontal dimensions. 

1. The Tank 

The tank, shown in Fig. 3.1, is made of wooden plates of 18 mm thickness 
which are sideways reinforced every 65 cm. The inner dimensions are 15.32 m 
length, 1.17 m width and 1.20 m depth. For better sealing the wood is covered 
with a layer of fiberglass. The tank rests on a concrete floor on foundation level. 

The Seawhisper™ anechoic tile system made by Seaward International, 
Inc. is positioned 8 cm away from the tank-walls. Two tiles are aligned vertically 
covering a height of the wall of 28 cm. The material reduces the echo reflection 
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by about 25 dB in the frequency range between 20 kHz and 100 kHz at normal 
incidence. An echo reduction versus frequency plot is shown in Fig. 3.2. The tank 
acts therefore as a horizontally infinite waveguide with almost rigid bottom and 
pressure release surface. 

On top of the tank sits a rail system on which two carts can move. The 
carts are able to carry equipment and hold the transducer arrays. One of them is 
driven by a step-motor, which is steered remotely from the computer system. 




Figure 3.1. View from above one end of the tank. The tank has inner dimension of 15.32 
m length, 1.17 m width and 1.20 m depth. The water level can be raised up 
to 28cm, the height of the anechoic material. 
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Lining System Performance: Rigid Backed 




Figure 3.2. Echo reduction versus frequency plot for the Seawhisper™ anechoic tile system 
covering the sidewalls. (A&T Engineering Technologies, VECTOR Research 
Division) 

2. The transducer arrays 

The two transducer arrays, made by EDO Electronic-Ceramic Products, 
consist of ten cylindrical piezo-ceramic elements each (see Fig. 3.3). The 
elements are positioned with its vertical cylindrical axis within a PVC tube. The 
dimensions of the arrays are 38 mm diameter and 305 mm length. 

The hoop mode frequency of all the elements is around 73 kHz. All 
elements by themselves in both arrays are horizontally omni-directional. The 
elements are all matched to one another very well as is shown in the calibration 
reports (EDO Corporation, Acoustic Division, September 30,1998). 
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305 mm 




Figure 33. The design of the ten-element transducer arrays made by EDO Electro-Ceramic 
Products of dimensions 38 mm diameter and 305 mm length. 



B. THE COMPUTERIZED DATA GENERATING AND 
ACQUISITION SYSTEM 

The data acquisition and generation system is made by Gage Applied 
Sciences Inc. and is PC-based. It consist of eight CompuGen 1100 function 
generator cards and four dual channel CompuScope 512 oscilloscope cards for a 
total of eight input and eight output channels available. 
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remote control 
transmitter 




preamplifiers 



computer 

system 



amplifiers 



Figure 3.4. The computer system with the preamplifiers, the 8-channel-switch, the amplifiers 
and remote control transmitter. 
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1. The computer and operating system 

The computer features a 266 MHz Pentium II CPU, 64 MB RAM, 6 GB 
hard-drive, 1 CD-ROM drive, 1 3.5” floppy disc, a VGA graphics card, a 15” 
monitor and a net card. The operating systems are DOS 7.0 and Windov\rs 98. 

2. The function generators 

The function generation is done with CompuGen 1100 function generator 
cards. Among their key features are: up to 80 Msamples per second, 12-bit 
resolution and up to 16 Msample memory buffers (Gage Applied Sciences Inc., 
February 1998). Function generator card number one acts as master and the 
seven other cards act as slaves. 

3. The oscilloscopes 

The data capture is done with CompuScope 512 oscilloscope cards. Their 
key features are: up to 5 Msamples per second, 12-bit resolution and 2 MB 
memory. They have all features of Digital Storage Oscilloscopes (Gage Applied 
Sciences Inc., February 1998). 

4. The software 

The cards are addressed via two MEX-files from Gage Applied Sciences, 
which are called by MATLAB programs. The MATLAB software was written by 
the author to fit our special needs. It includes the key features definition of 

signals to be transmitted, editing of captured signals, including cutting, amplifying 
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and time-reversing, setting the parameters of and initializing the data system, 
data storing, visualization and analysis of the data. The software offers four 
graphical user interface windows, as shown in Fig. 3.5. The main window allows 
the experiment to run, while the two signal creation windows allow prepare 
signals to be prepared and messages to be used in the experiment. The remote 
control window allows changing the position of one equipment cart and both 
arrays in order to change the geometry of the experiment. 

5. Amplifier 

To amplify the output of the function generators we used a set of eight 
HP 467A amplifiers. This allowed us to amplify the signals up to ±20 V peak to 
peak. 



6. Preamplifier 

Pre-amplification of the received analog signals, before they were 
captured by the oscilloscope cards, was done with a set of eight Stanford SR560 
preamplifiers. Preset 6dB/ocatve filtering of the incoming signals can be 
accomplished directly with the built-in filters. 

7. Switch 

To switch the eight outgoing and the eight incoming channels from a 
parallel (meaning amplifiers to transducer#! and transducer #2 to preamplifiers) 
to a crosswise (meaning amplifiers to transducer #2 and transducer#! to 
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preamplifiers) connection we used a relays switch arrangement, using a total of 8 
relays (one per channel). 

8. Remote Control Equipment 

To remotely move the position of one equipment cart in range and the two 
transducer arrays in height, we used a radio transmitter set and stepping motors. 
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IV. EXPERIMENT DESCRIPTION 



In this chapter we describe the experiments we conducted using TRA 
techniques as applied to underwater acoustic communications. 

Using common procedures for calibration and testing, we found it best to 
use a frequency range between 50 kHz and 65 kHz. Below this range, the 
response of our system was weak, and electronic noise at 30 kHz and 35 kHz 
would yield a small signal-to-noise ratio (SNR); also switching the relays into the 
on position would cause 60 Hz noise. Above the frequency range of operation, 
the strong response at resonance (~73 kHz) overshadowed the response at 
other frequencies. 

A. DEFINITION OF SYMBOLS 

In the following we will use the expression symbol for a single pulse 
containing one, two, three or all four center-frequencies. A symbol thus contains 
information which can be related to a code. Four center-frequencies can be 
combined into 15 different ways, excluding the null symbol which contains no 
acoustic energy. For each combination we assigned a number between 1 and 
1 5 as shown in Table 4. 1 . 



SYMBOL NUMBER 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


13 


14 


15 


Center-Frequency 1 


X 




X 




X 




X 




X 




X 




X 




X 


Center-Frequency 2 




X 


X 






X 


X 






X 


X 






X 


X 


Center-Frequency 3 








X 


X 


X 


X 










X 


X 


X 


X 


Center-Frequency 4 
















X 


X 


X 


X 


X 


X 


X 


X 



Table 4.1. The center-frequencies used in each symbol. 
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B. RESOLUTION OF SYMBOLS 



The four center frequencies we used are at 50 kHz, 55 kHz, 60 kHz and 
65 kHz. The initial pulse transmissions used one single center frequency at a 
time. The pulses were created using different bandwidths, corresponding to 
different lengths in time. The pulses were then multiplied by a Hanning window, 
but we also used a rectangular and a Bartlett window during our testing. 

1. The initial transmission 

The initial transmission commonly was done in the parallel switch setting, 
meaning that transducer array #1 transmitted and transducer array #2 received. 
The resulting reception contained all the multipath propagation. We then selected 
a window that contained all the multipath propagation information to create the 
time reversed signal, making sure that the window size within one set of data 
was the same and that the relative position of the time windows at the first arrival 
was the same. 

2. The messages 

The operation described above resulted in time reversed signals which we 
used to build the messages containing various symbols. Usually we started with 
a message using the symbols ‘1-2-4-8-15’. The spacing between the symbols 
was varied. This allowed discrimination of the response of the single frequencies 
as well as the response of the combination of all of them. For further messages 
we could use the knowledge about the response to adjust the amplitudes 
between the center-frequencies relative to each other. The message we defined 
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in this way was sent with the crosswise switch setting, meaning that transducer 
array #2 transmitted and transducer array #1 received. 

3. Bandwidth limit and spacing between symbols 

We send the messages created as described above and checked whether 
we could resolve them by a matched filter approach and by a FFT approach. We 
did this for bandwidths of 10 kHz, 5 kHz, 4 kHz, 3 kHz, 2.5 kHz, 2 kHz, 1.5 kHz, 
and 1 kHz. The spacing between symbols was also varied from none to one 
symbol length. 



C. FOCAL REGION SIZE 
1. Full array 

To determine the focal region size we used a single time-reversed signal 
rather than a symbol or a message. Starting at the range at which the time- 
reversed signal was received we decreased and increased the range by steps of 
1 cm. The time-reversed signal was sent at each of the positions and the signal 
received was checked for time and spatial (in depth) focus. We 
increased/decreased the range until time-focus was not resolved. 
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2. Reduced number of elements 

For an amplitude adjusted time-reversed signal with frequency 50 kHz and 
bandwidth 2000 Hz we used one, two and three element combinations to 
transmit the time reversed signal. We compared these receptions with the all- 
element array reception for the same signal. This comparison allowed us to 
determine the minimum number of elements that would still give good temporal 
and spatial focus. 



D. FOCAL REGION STEERING IN RANGE 

We shifted the frequencies of the time-reversed signals in order to focus at 
different ranges. These new time-reversed signals were transmitted as before, 
with the receiving transducer moved to the new expected focus location. Then 
we increased/decreased the range in steps of a wavelength in order to find the 
size of the new focal region and its center. The expected focal range was 
determined from Eq. (2.24), with a value oi P = 1.0 . This turned out to be a good 
starting point. 
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V. RESULTS 



In this chapter we will describe and analyze the results of the experiments 
described in chapter IV. 

A. RESOLUTION OF SYMBOLS 

The resolution of symbols depends on the quality of the reception and the 
means used to resolve the frequency data they contain. We used two simple 
approaches. In the first we performed an FFT on each symbol separately. In the 
second approach we analyzed the symbols by match filtering the received 
signals with the individual initial signals at each center frequency. 

1. Bandwidth and spacing limits using a FFT 

The three quantities that determine data rates are inter-symbol spacing, 
bandwidth, and center frequency spacing. The bandwidth determines the time 
span of the symbol, and the frequency spacing determines the resolution of the 
information contained in a symbol. The inter-symbol spacing is mainly 
determined by the energy of the sidelobes at the focal point. The use of a 
Hanning window minimized the sidelobes allowing also for no spacing between 
the symbols. This is not the case for a rectangular window. Thus, under the 
condition of no symbol spacing, it is important to select frequency signals which 
give a reception that can be resolved in its frequency components. This condition 
was fulfilled for the center frequencies mentioned previously with a bandwidth of 
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2 kHz. Fi ire 5.1 shows a message with 2.5 kHz bandwidth, which is not 
resolvable . id one with 2 kHz bandwidth, which is clearly resolvable. Note that 
because the center frequency spacing is 5 kHz, there exists a ratio of 5 to 2 for 
the center frequency spacing to the bandwidth. 

We set the thresholds for accepting the signal after undergoing a FFT as a 
minimum of 6 dB difference between the lowest amplitude valid peak and the 
highest noise peak. Also we required that the valid peaks all lay within 9 dB of 
the highest amplitude within the symbol. We found these conditions necessary in 
order to have automation capability for analyzing the messages. 

We also varied the spacing between symbols. The spacing was 
decreased from a full symbol length to no spacing at all. In Fig. 5.2a and Fig. 
5.2b we show the FFT for a spacing of one symbol length and no spacing 
respectively . 

In the rest of the analysis we will only use messages without symbol 
spacing, since this is the more desirable case in order to have high data rates. 
Also we will use Hanning window signals throughout. 
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signal level [db] 



Bandwidth 2500 Hz 




Frequency [kHz] 
Bandwidth 2000 Hz 







Frequency [kHz] 

Figure 5.1 Two messages containing the symbols 1 - 2 - 4 - 8. The left message has a 

bandwidth of 2500 Hz and is not resolvable. The right message has a bandwidth 
of 2000 Hz and it is clearly resolvable. 
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Signal level [db] 







Frequency [kHz] 




Frequency [kHz] 





Figure 5.2.a Two messages with one symbol length spacing. The left message contains the 
symbols 1 - 2 - 4 - 8 and the right message the symbols 15-9-6. 
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Signal level [db] 






Frequency [kHz] 




Frequency [kHz] 





Figure 5.2.b Two messages with no spacing. The left message contains the symbols 1-2-4 
- 8 and the right message the symbols 15-9-6. 
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2. Bandwidth limits using a Matched Filter 

For a message containing Hanning window signals and having no spacing 
we found that there is a maximum bandwidth of 2500 Hz, when the received 
symbols were match filtered with the initial signals whose frequency spacing is 5 
kHz. Therefore we find a ratio of 2 to 1 for center frequency spacing to 
bandwidth. 

Our criteria for accepting a signal as resolvable is that the sidelobes must 
be less than half the amplitude of the main peaks. This criteria is based on the 
notion that when the sidelobes of two adjacent symbols are added, the symbol 
can still be resolved. 

Figure 5.3 shows messages with 2.5 kHz and 3 kHz bandwidth. The 2.5 
kHz bandwidth message barely meets our criteria, while the 3 kHz bandwidth 
messages does not. 

The matched filter approach to resolve the symbols of the messages 
turned out to be the more desirable one, both because the resolvable bandwidth 
is larger and also because it is easier to apply. With the matched filter the 
frequency content of all symbols in a message is deciphered at once without 
having to define the time window for execution as it is done with the FFT 
approach. 
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Matched Filtering Result - Bandwidth 2500 Hz 



50000 Hz 



55000 Hz 



60000 Hz 



65000 Hz 




3 4 

Time [msec] 



Matched Filtering Result • Bandwidth 3000 Hz 




0 0.5 1 1.5 2 2.5 3 




Figure 5.3. Two messages containing the symbols 1 - 2 - 4 - 8. The left message with a 

bandwidth of 2500 Hz is resolvable, the one on the right with 3000 Hz bandwidth 
is not resolvable. 
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B. FOCAL REGION SIZE 



1. Full array 

We used 2kHz bandwidth signals centered at 50 kHz and 65 kHz to 
determine the extent of the focal region. We found that the focal region has an 
ellipsoidal cross-section, with tt>e major axis along the range direction and the 
minor axis along depth. 

The extent in depth of the focal region is ± 15 cm as measured at the half- 
power point. This corresponds to a distance of about 0.6i . The simulations for a 
frequency range of 400 Hz to 1200 Hz showed about focal region extend. 
Figure 5.4 shows the normalized, curve-fitted amplitude distribution over the 
depth for a signal that was sent from element two. The depth is shown relative to 
the intended focal position. 

It can be seen that the energy distributions for both frequencies differ. The 
higher frequency has a stronger sidelobe. This is due to a higher attenuation of 
the higher modes contained in the higher frequency signal, which causes the 
sidelobes to become stronger. This is in contrast to the simulation results done at 
frequencies of 400 Hz to 1200 Hz, where the higher frequencies had the smaller 
focal region size. The focal point is shifted slightly for the two frequencies, where 
the higher frequency is shifted a bit downward and vice versa. Both frequencies 
show an increased amplitude at depths close to the rigid bottom. This is a 
phenomena to be expected near a rigid bottom and has no relation to the focal 
region. 
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- 50 kHz 
-65 kHz 




1 ' 1 ' 1 ' 1 ' 1 ' 1 ^ 1 ' 1 

-2 0 2 4 6 8 10 12 



Relative Depth [cm] 



Figure 5.4. Focal region size in depth for a 50 kHz and a 65 kHz signal. The depth is 
measured relative to the intended focal point. 

The extent in range of the focal region is much larger. For 50 kHz we 
found +9.5 cm and -8 cm at the half-power point. This is equivalent to about ± 3 
wavelengths. The focal region for the 65 kHz signal extends +1 1 cm in direction 
of increasing range, slightly more than the 50 kHz signal. However, it is shifted in 
decreasing range by -5 cm and extends to -20 cm in the decreasing range 
direction. The equivalent wavelengths are about +5A and -9A. Again the 
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simulations for a frequency range of 400 Hz to 1200 Hz showed a significantly 
larger extend. 

The shifting and the sidelobe phenomena already noticed for the focal 
region along depth becomes more pronounced along range. The loss of energy 
and therefore information in the higher modes contained in the higher 
frequencies leads to a degradation of the focal region by a significant amount. 
The results are shown in Fig. 5.5. 




Figure 5.5. Focal region size in range for a 50 kHz and a 65 kHz signal. The range is 
measured relative to the intended focal point. 
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2. Focal region for 1 element 



The response for one channel transmissions showed a weak focusing in 
space for most transmission channels as shown in Fig. 5.6. However, the time 
focus was still good for the intended focal point. It therefore can be concluded 



Full Array 




1 Channel Series 




Relative Depth [cm] 



Figure 5.6. Focal region size in depth for a 50 kHz. The left side shows the full array 

response, the right side the response for one transmission channel for each of 
the eight channels. The depth is measured relative to the intended focal point. 
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Channel 



that, if only a time focus is wanted, it is good enough to use a single element. 
The response will then be strongly dependent on location of the transmission. 
The closer the transmitting transducer is to a layer boundary, the surface or the 
bottom, the worse the response will be. 




Time [ms] 



Figure 5.7. Time focus versus depth for a 50 kHz signal send only from channel 3. The 

energy is spread both in time and depth. However, for the intended focal point on 
channel five the time focusing is acceptable. 
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Channel 




0 0.20002 0.40004 0.60006 0.80008 1.0001 1.20012 1.40014 1.60016 1.80018 2.0002 



Time [ms] 

Figure 5.8. Time focus versus depth for a 50 kHz signal send from all channels. The energy 
is well focused on the intended focal point. 

In Fig. 5.7 the intensity in time versus depth (channels) of a transmission 
from channel 3 is shown. This is to be compared to the response of the full array 
transmission as shown in Fig. 5.8. 
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Normalized [Amplitudel Normalized |Amplitude| 



3. Focal region for 2 elements 



The response for two channel transmissions showed a significantly better 
focusing in space for most transmission channels, as shown in Fig. 5.9, than 



Full Array 




2 Channel Series 




Relative Depth [cm] 



Figure 5.9, Focal region size in depth for a 50 kHz. The left side shows the full array 

response, the right side the response for two successive transmission channels 
for each of the seven possible combinations. The depth is measured relative to 
the intended focal point. 
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Channel 



the 1 channel transmissions. The time focus was also better for the intended 



focal point as is shown in Fig. 5.10 (again compare to Fig. 5.8). 




Time [ms] 



Figure 5.10, Time focus versus depth for a 50 kHz signal send from channel 3 and 4. The 

energy is spread both in time and depth. However, there is a good focus in time 
given for the intended focal point on channel five. 
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4. Focal region for 3 elements 

The response for a three channel transmission showed a very good 
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Figure 5.11. Focal region size in depth for a 50 kHz. The left side shows the full array 

response, the right side the response for three successive transmission channels 
for each of six possible combinations. The depth is measured relative to the 
intended focal point. 
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Channel 



focusing for both time and space as can be seen in Fig. 5.11 and Fig. 5.12. This 
should be expected. The more elements of a vertical array used, the better will 
be the focusing. Thereby three array elements within a wavelength are 
commonly used to get good sampling. However, an array with two elements, 
arranged with a distance of less than a wavelength is sufficient, as we could 
show with our data. 



Intensity Plot Of 3 Channel Transmission - Channels 3-5 




Time [ms] 



Figure 5.12. Time focus versus depth for a 50 kHz signal send from channel 3,4 and 5. The 
energy is spread in depth. There is a very good focus in time for the intended 
focal point on channel five and a satisfying time focus on the other channels. 
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C. FOCAL REGION STEERING IN RANGE 



By shifting the frequency, the focal region shifts. The measurements for 
the shifted focal region showed that the simple waveguide we used has a value 
of y? = l. For the 50 kHz signal a shift in range of 20 cm is achieved by a 
frequency shift of 2500 Hz. The focal region extent was the same for both the 
unshifted and shifted signal as is shown in Fig. 5.13 and Fig. 5.14. We conclude 
that shifting the focal region by 5% in range does not significantly change the 
size of the focal region. 




Figure 5.13. Focal region size in range of the unshifted and shifted signal. Both are measured 
and shown relative to the intended focal point. 



48 




Figure 5.14. Focal region size in depth of the unshifted and shifted signal. Both are measured 
and shown relative to the intended focal point. 



D. MESSAGES 

Because the focal region is limited in size and there are strong sidelobes 
outside of that region, there is a natural encryption built into the transmission of 
messages. The signals do not focus well in time outside of the focal region and 
therefore there is inter-symbol interference that makes it impossible to recover 
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the messages. In order to prove this we defined a simple code using the 15 
symbols we get from the 4 center-frequencies. This code is shown in Table 5.1. 



General 


Symbol 


Meaning 


14 


number follows 


15 


word follows 


_ (space) 


shift for second half of letters 



Letters 


Symbol 


Meaning 


1 


A 


2 


B 


3 


C 


4 


D 


5 


E 


6 


F 


7 


G 


8 


H 


9 


1 


10 


J 


11 


K 


12 


L 


13 


M 




1 


N 


2 


0 


3 


P 


4 


Q 


_5 


R 


6 


S 


7 


T 


8 


U 


9 


V 


10 


w 


11 


X 


12 


Y 


13 


z 



Numbers 


Symbol 


Meaning 


1 


1 


2 


2 


3 


3 


4 


4 


5 


5 


6 


6 


7 


7 


8 


8 


9 


9 


10 


0 


11 


+ 


12 


- 


13 





Table 5.1. Simple code using 15 symbols. 

1. Message at the intended reception point 

The message we used was “DEC 99”, which using the simple code above 
reads “15-4-5-3-14-9-9” in symbols. This message could clearly be read at the 
focal point using a matched filter analysis. Figure 5.14 shows the matched filter 
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result. Comparing this with the Table 5.2, in which the center-frequencies of the 
symbols transmitted are shown, shows that the message is correctly transmitted 
and easily resolved. 



15 4 5 3 14 9 9 




0 1 2 3 4 5 6 

Time [msec] 



Figure 5.14. Matched-filter analysis of the incoming message. The message is clearly and 
easily resolvable. The red boxes show the symbols with respect to time. 



Sym bol 


1 5 


4 


5 


3 


14 


9 


9 


50 kHz 


X 




X 


X 




X 


X 


55 kHz 


X 






X 


X 






60 kHz 


X 


X 


X 




X 






65 kHz 


X 








X 


X 


X 



Table 5.2. Center-frequencies transmitted for the symbols of the message “DEC 99” using 
the simple code defined in Table 5.1. 
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2. Message outside the focal region 

The message outside the focal region is naturally encrypted. This can be 
seen in Fig. 5.15, where the message from Table 5.2 is used again but is 
received outside the focal region. 



15 4 5 3 14 9 9 




0 1 2 3 4 5 6 



Time [msec] 



Figure 5.14. Matched-filter analysis of the incoming message outside the focal region. The 
message is unclear and not resolvable. The red boxes show the symbols with 
respect to time. 



Svm bol 


14 
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2 


1 4 


1 0 


8 


50 kHz 
















55 kHz 


X 






X 


X 


X 




60 kHz 


X 


X 


X 




X 






65 kHz 


X 








X 


X 


X 



Table 5.3. Center-frequencies as received outside the focal region resulting in the symbols 
of a message reading “442 08” using the simple code defined in Table 5.1 . 
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Note that for the higher frequencies the signal becomes clearer and easier 
to resolve. This is due to the higher attenuation of the higher modes of the partial 
pulse. The energy propagating on the direct path is less attenuated and 
therefore relative to the other paths is dominating. This reduces the effect of 
natural encryption. Also, it should be noted that the response of our system is 
less strong for the lower frequencies due to the frequency roll-off of the elements 
of the array. 

For the message shown we used a bandwidth of 2000 Hz. This is 
equivalent to a data rate of 2000 symbols per second. Since the symbols we 
used where built from four center-frequencies we can also express that number 
as 8000 bits per second. We have shown earlier that the limit we could achieve 
with our system was a bandwidth of 2500 Hz. This is equivalent to a limit of 2500 
symbols per second or 10000 bits per second. 
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VI. CONCLUSIONS 



A. CONCLUSIONS 

In this thesis experimental studies were conducted to test results of 
analytical and numerical work on time-reversed signal underwater 
communication. The experimental results indicated a good basis for further 
experimental work both in scale model experiments and at-sea experiments. 

The experimental results agree with earlier numerical work in regards to 
the focusing properties. That is, the size of the focal region is larger in range 
(horizontal) than it is in depth (vertical). The focal region size was found to be 
only a few wavelengths along the horizontal direction and less than one 
wavelength along the vertical direction. Outside the focal region intersymbol 
interference occurs due to the lack of time focusing, which provides a natural 
encryption property. This can be of advantage especially in military use. It also 
means that a vertical change of position is more critical to communication than a 
horizontal change. 

We have also shown that the focal region can easily be shifted in range by 
frequency shifting of the time-reversed signal. If a planar array is used, a shift in 
the bearing angle can be achieved with commonly used beamforming 
techniques. Thus, the location of the focal region can be adjusted to stay at the 
position of a moving receiver, such as an autonomous underwater vehicle or 
submarine. 
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Experiential results also show that there is an optimal value for the relation 
between bandwidth and spacing of center-frequencies. For 5 kHz center- 
frequency spacing the optimal bandwidth was found to be to be 2.5 kHz. Thus, a 
total bandwidth of 10000 Hz is needed to transmit a data rate of 10000 bits per 
second. These results show that the TRA approach offers powerful means for 
high data rate underwater communication because of reduced inter-symbol 
interference and bandwidth efficiency. 

A problem arising from the stronger attenuation of higher modes of the 
time-reversed signal is the degradation in spatial focusing. The time focusing was 
however not affected from this attenuation. 

The spatial focusing properties depend on the aperture size and number 
of elements per wavelength. This work showed a slight improvement in the 
spatial side-lobe amplitude for the use of a larger aperture while keeping the 
element spacing constant. The degradation of the focal region with less aperture 
size and therefore less number of elements is due to reduced spatial sampling as 
well as less energy transmitted. However, focusing in time still occurred, 
indicating that the necessity of guard times with a time-reversal approach can be 
eliminated. 



B. RECOMMENDATIONS FOR FURTHER EXPERIMENTS 

Steering the focal region could only be investigated for a range shift. 
However, a shift in relative depth needs to be investigated. In general, depth 
shifting of the focal region can be done by shifting the relative phase between the 
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elements of the array. Note that the phase shift has to be applied to a signal 
containing multimode information, making this a very involved and difficult 
process. 

Further efforts should be made in investigating the application of other 
communication approaches in order to combine them and thereby increase the 
data rate. For example, the use of different phase information for a fixed center 
frequency might be a way to combine the algorithm and increase the data rate. 

Before this approach can be implemented in a real system, it is necessary 
to prove the reliability in real environments. A step in that direction could be to 
add noise in a simple environment like our tank and to increase the complexity of 
our environment. 

Finally, a very interesting and instructive area of research is to perform 
benchmark testing of different communication schemes and determine their 
limitations and advantages over the TRA technique. A group of researchers at 
the Naval Postgraduate School is already embarking on such a task. 
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